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Orthostatic Intolerance (OI) 


Propensity to develop symptoms of fainting during 
upright standing. 

OI is associated with changes in heart rate, blood 
pressure and other measures of cardiac function. 



Orthostatic Intolerance (OI) 


Propensity to develop symptoms of fainting during 
upright standing. 

OI is associated with changes in heart rate, blood 
pressure and other measures of cardiac function. 


Problem: NASA astronauts have shown increased 
susceptibility to OI on return from space missions. 


Current methods for counteracting OI in astronauts include 
fluid loading and the use of compression garments. 
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Assessment of OI: Orthostatic Tilt Tests (OTTs) 

• Subject initially is supine. 

• 80° upright tilt for preset time 
(T max = 5-30 min.) 

• “Survival” time = T 

• Endpoint: T c = min(7; T max ) 



Concurrent measurements x(t ) = [x ; (/),..,x s (t)] (t <= T c ): 

Xj{t) = heart rate (hr) 

x 2 (t)...x 5 (t) = measures of blood pressure ( dbp , map, sbp, pp) 

x 6 {t) = stroke volume (.vv) 

x 7 (t) = cardiac output ( co = sv x hr) 

x 8 (t) = total peripheral resistance tpr = (map - mvp)/co 



100 50 100 150 200 


Time trajectories ofx It ..,x 8 (completed OTT) 



hr — • — dbp • sbp • — map • pp 



Xj heart rate (hr) 

x 2 disastolic bp (dbp) 

x 3 mean arterial bp (map) 

x 4 systolic bp (sbp) 

x 5 pulse pressure (pp) = sbp-pp 

x 6 stroke volume (sv) 

x 7 cardiac output (co) 

x 8 total peripheral resistance (tpr) 
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Uncompleted OTT 

Time trajectories of [x/t); t = 1 ,..,T C } 



Causative Flow (no censoring) 

We assume all the information about the degree of 01 would be 
contained in the survival time T, if there were no censoring. 


01 — *■ Concurrent measurements {x(t);(t<T)} — > T 
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Causative Flow (censoring) 

With censoring present, both (t < T c )} and T c provide 
information about the survival time T. 

WO; (!<T C )} 

T c = min ( T, T max ) 
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Bedrest simulation of spaceflight 


• 27 subjects (10F, 17M) 

• 60-day bedrest 

• OTTs pre and post-bedrest 

• r max = 30min - 

Endpoint: T c = min( T, T max ) (in min.) 

concurrent measurements : (t = 0, 1,2, ... ,T C ) 

Xj(t ) = heart rate (hr) 

x 2 (t)...x 5 (t) = measures of blood pressure ( dbp , map, sbp, pp) 

x 6 (t ) = stroke volume (sv) 

x 7 (t) = cardiac output ( co = sv xhr) 

x 8 (t) = total peripheral resistance tpr = (map - mvp)/co 
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OTT survival times for bedrest subjects, by gender 


Female 


id 

pre BR 

post BR 

4 
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8 

30 

20 

9 
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11 

12 
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7 
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6 

8 

21 

10 


22 

22 



id 

pre BR 

post BR 

1 

30 

30 

2 

18] 

8 

3 

30 

22 

5 

28 

24 

6 

20 


7 

30 


12 

30 


13 

28 


14 

30 


16 

22 

11 

19 

19 

29 

20 

30 

3 

23 

17 

6 

24 

30 

1 

25 

30 

9 

26 

12 

3 

27 

30 

26 


Male 
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Pre vs. post bedrest comparison (univariate) 


Sub 1 


Pre BR Post BR 



Graphs by brday 


Sub 2 


Pre BR Post BR 



Xj = heart rate (hr) 

x 2 = diastolic blood pressure (dbp) 

x 3 = mean arterial pressure (map) 


Sub 3 


Pre BR Post BR 



Graphs by brday 


Graphs by brday 
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2-d Trajectories 


raw values with supine time point 


Sub 1 



DBP 
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2-d Trajectories 


raw values with supine time point 

Sub 2 

Pre BR Post BR 


o 
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Graphs by brday 


dbp 



DBP 


100 150 200 


2-d Trajectories 


raw values with supine time point 


Sub 3 



DBP 


Summarizing Behavior of Concurrent Measurements 

In practice we seek a summary statistic A that captures the 
essential information in (x(t)}. 

A=A(x(t)) ( t<T c ) 

y 

oi 

^ T C = min ( T, T max ) 


Desirable properties of A 

• Summarizes relevant behavior of (X](t),..jc 8 m 

• Can be calculated for OTT of any (reasonable) length. 

• A and T c explain 01 better than T c alone for short tests. 
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Candidates for A = A(x(t)) 

/^-dimensional functional data x(t) is observed at discrete common 
time points 


1 . determinant: A k = \S k \ ]/ p ; S k 



k 


ll 


-x k )(x(ti) -x k y 


k 

2. path length: ^=^1^*0 -^*(h-i)l 

L = 1 


3. convex hull area/volume: A k 
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Properties of A=A(x(t)) 


Generally increases with time. 


isub=3 

Pre BR Post BR 



0 10 20 30 0 10 20 30 

minute 


Graphs by brday 


For longer surviving subjects, average slope is lower. 
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Properties of A=A(x(t)) 


Generally increases with time. 


isub=3 

Pre BR Post BR 



0 10 20 30 0 10 20 30 

minute 


Graphs by brday 


For longer surviving subjects, average slope is lower. 



Outcome for inference: y = A(x(TJ)/T c 
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Inference on effect of bed-rest on OI. 


Vy — // + Uj + OCj + €■■ (7-th subject; /-th treatment) 


T < T 

ij — max 

T- > T 

ij max 


}'ij = U.V( /,))//, 

y tj e ( 0, A(x(T max ))/T max ) 


CO - 

T- = survival time (may be censored) 


CM - 



0 


K) 

tmax 


Assume ~ Normal 

max likelihood (integrate out random effects) 
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Exploratory Study Plan 


Observe survival and functional data x(t ) from OTTs 
given to subjects pre- and post-bedrest. 

Formulate some candidates for A = A(x(t)) (t < T c )). 

Compare their ability to test for an effect of bedrest at 
various censoring times with a simple non-parametric 
analysis. 
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Cox regression (OTT survival times) 

censoring time: minutes 


. stcox brday, cluster (isub) nolog 

failure _d: fail 

analysis time t: tmax 


Cox regression -- Breslow method for ties 


No. 

of subjects = 

46 

Number 

of obs = 

46 

No. 

of failures = 

34 




Time at risk = 

840 







Wald chi2 (1) 

17.03 

Log 

pseudolikelihood = 

108.11632 

Prob > 

chi2 = 

0.0000 



(Std. 

Err. adjusted for 

27 clusters 

in isub) 


1 

Robust 





t | Haz . Ratio 
1 

Std. Err. 

z P> | z | 

[95% Conf. 

Interval] 


T 

brday | 3.189786 

.8966697 

4.13 0.000 

1.838583 

5.534009 
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Cox regression (OTT survival times) 

censoring time: minutes 


. stcox brday, cluster (isub) nolog 

failure _d: faill5 

analysis time _t: tmaxl5 

Cox regression -- Breslow method for ties 

46 

17.29 
0.0000 


No. of subjects 
No. of failures 
Time at risk 


46 

19 

551 


Log pseudolikelihood = -63.369599 


Number of obs 


Wald chi2 (1) 
Prob > chi2 


(Std. Err. adjusted for 27 clusters in isub) 


1 

t | Haz . Ratio 
1 

Robust 
Std. Err. 

z 

P>|z| 

[95% Conf. Interval] 

T 

brday | 

1.730109 

4.16 

0.000 

2.258201 9.652309 
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Cox regression (OTT survival times) 

censoring time: 10 minutes 


. stcox brday, cluster (isub) nolog 

failure _d: faillO 

analysis time t: tmaxlO 


Cox regression -- Breslow method for ties 


No. 

of subjects = 

46 

Number 

of obs = 

46 

No. 

of failures = 

13 




Time at risk = 

406 







Wald chi2 (1) 

13.00 

Log 

pseudolikelihood = 

■43.233908 

Prob > 

chi2 = 

0.0003 



(Std. 

Err. adjusted for 

27 clusters 

in isub) 


1 

Robust 





t | Haz . Ratio 
1 

Std. Err. 

z P> | z | 

[95% Conf. 

Interval] 


i 

brday | 6.313793 

3.227091 

3.61 0.000 

2.318602 

17.19311 
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Cox regression (OTT survival times) 

censoring time: 7 minutes 


. stcox brday, cluster (isub) nolog 

failure _d: fail7 

analysis time _t: tmax7 

Cox regression -- Breslow method for ties 

46 

4.98 
0.0257 


No. of subjects 
No. of failures 
Time at risk 


46 

8 

300 


Log pseudolikelihood = -25.673825 


Number of obs 


Wald chi2 (1) 
Prob > chi2 


(Std. Err. adjusted for 27 clusters in isub) 


1 

Robust 





t | Haz . Ratio 
1 

Std. Err. 

z 

P>|z| 

[95% Conf. 

Interval] 

i 

brday | 11.96172 

13.3078 

2.23 

0.026 

1.351463 

105.8725 
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Cox regression (OTT survival times) 

censoring time: minutes 


. stcox brday, cluster (isub) nolog 

failure _d: fail5 

analysis time _t: tmax5 

Cox regression -- Breslow method for ties 

No. of subjects = 46 Number of obs = 46 

No. of failures = 

Time at risk = 221 


[no analysis] 
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Best Results 

(by Outcome Type and Censoring Time) 


cens time 

path length 

Z 

convex hull 

Z 

determinant 

Z 

5 

HR*, DBP* 

4.97 

HR 

4.71 

HR, MAP, SV 

3.05 

7 

HR*, PP* 

6.59 

HR, DBP, SBP 

4.48 

HR, DBP, SBP 

4.5 

10 

HR*, TPR*,PP* 

7.13 

HR, DBP, SBP 

6.08 

HR, DBP, SBP 

5.84 

15 

HR*, TPR*,PP* 

7.58 

HR, DBP, SBP 

6.67 

HR, DBP, SBP 

6.05 

30 

HR*, TPR*,PP* 

6.65 

HR, DBP, SBP 

6.49 

HR, DBP, SBP 

6.28 
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Summary of approach used to developing 
the statistic A = A(x(t )); t < T c 


1 . Used training data with control condition and 
provocative condition known to cause 01. 
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Summary of approach used to developing 
the statistic A = A(x(tJ); t < T c 

1 . Used training data with control condition and 
provocative condition known to cause 01. 

2. Observed increased variability of functional 
data x(t ) for the provocative condition. 
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Summary of approach used to developing 
the statistic A = A(x(tJ); t < T c 

1 . Used training data with control condition and 
provocative condition known to cause 01. 

2. Observed increased variability of functional 
data x(t ) for the provocative condition. 

3. Formulated several measures of this variability 
as candidates for A. 
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Summary of approach used to developing 
the statistic A = A(x(t )); t < T c 

1 . Used training data with control condition and 
provocative condition known to cause 01. 

2. Observed increased variability of functional 
data x(t) for the provocative condition. 

3. Formulated several measures of this variability 
as candidates for A. 

4. Used simple non-parametric test that combines 
survival data with cumulative evaluations of 
each ^-candidate to make inference on the 
effect of condition. 
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Summary of approach used to developing 
the statistic A = A(x(t )); t < T c 

1 . Used training data with control condition and 
provocative condition known to cause 01. 

2. Observed increased variability of functional 
data x(t) for the provocative condition. 

3. Formulated several measures of this variability 
as candidates for A. 

4. Used simple non-parametric test that combines 
survival data with cumulative evaluations of 
each ^-candidate to make inference on the 
effect of condition. 

5. Compared results. 



conclusions / remarks 


Multivariate trajectory spread is greater as 01 increases. 

Pairwise comparisons at the same time within subjects 
allows incorporation of pass/fail outcomes. 

Path length, convex hull area, and covariance matrix 
determinant do well as statistics to summarize this spread 

Missing data problems 

Time series analysis 

need many more time points per OTT session 
treatment of trend? 

how incorporate survival information? 
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